Evidence of silicene in honeycomb structures of silicon on Ag(lll) 
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In the search for evidence of sihcene, a two-dimensional honeycomb lattice of silicon, it is important 
to obtain a complete picture for the evolution of Si structures on Ag(lll), which is believed to be 
the most suitable substrate for growth of silicene so far. In this work we report the finding and 
evolution of several monolayer superstructures of silicon on Ag(lll) depending on the coverage and 
temperature. Combined with first-principles calculations, the detailed structures of these phases 
have been illuminated. These structure were found to share common building blocks of silicon 
rings, and they evolve from a fragment of silicene to a complete monolayer silicene and multilayer 
silicene. Our results elucidate how silicene formes on Ag(lll) surface and provide methods to 
synthesize high-quality and large-scale silicene. 
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Introduction With the development of semiconduc- 
tor industry toward smaller scale, the rich quantum phe- 
nomena in low-dimensional systems may lead to new 
concepts and ground-breaking applications. In the last 
decade graphene has emerged as a low-dimensional sys- 
tem for both fundamental research and novel applications 
including electronic devices, energy storage and transpar- 
ent protection layer [l-|j| . Inspired by the fruitful results 
based on graphene, recently a lot of interest has been 
drawn to group IV (Si, Ge) analogs of graphene. [5|- 
[tI] It has been theoretically shown that silicene, with Si 
atoms packed in a honeycomb lattice like graphene, is 
a new massless Dirac Fermion system. [5|, |8| Compared 
with that of graphene, the stronger spin-orbit coupling 
in silicene may lead to detectable quantum spin Hall ef- 
fect (QSHE) and other attractive properties. [8|-[ll| The 



compatibility of silicene with silicon-based nanotechnol- 
ogy makes this material particularly interesting for device 
applications. 

As the theoretical studies on silicene are rapidly in- 
creasing, the major challenge in this field is now the 
preparation of high quality silicene films. However, to 
date, there is still no solid evidence for the observation of 
a silicene film. There have been a few works on the forma- 
tion of silicene nanoribbons on Ag(llO) with graphene- 
like electronic signature 12|, [l3|. The only published 



work on the preparation of silicene-like sheets was re- 
ported by Lalmi et al. on Ag(lll) [IJj. They showed 
scanning tunneling microscopy (STM) images of hon- 
eycomb monolayer structure that resembles monolayer 
graphene structure. However, in their experiment the 
observed lattice constant was about 17% smaller than 
the theoretically proposed model or the value of bulk sil- 
icon. Such a huge compression of the lattice is rather 
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unlikely to be induced by the strain between the film 
and the substrate. Their results therefore remains to be 
confirmed and understood. For the purpose of finding 
evidence of silicene, and optimizing the preparation pro- 
cedure for growing high quality silicene films, it is impor- 
tant to build a complete understanding of the formation 
mechanism and growth dynamics of possible silicon struc- 
tures on Ag(lll), which is currently believed to be the 
best substrate for growing silicene. 

In this paper, we present a systematic study of the 
self-organized superstructures formed by sub-monolayer 
silicon grown on Ag(lll), by STM and scanning tunnel- 
ing spectroscopy (STS). We found that, depending on 
the substrate temperature and silicon coverage, several 
monolayer superstructures can form on Ag(lll). These 
superstructures are distinct from any known surface 
structures of bulk silicon, and are characterized by honey- 
comb building blocks and structures. At sufficiently high 
temperature and Si coverage, monolayer and multilayer 
silicene films were grown. Combined with first-principles 
calculations, the structural models of these phases are 
proposed and their evolution with temperature and Si 
coverage is discussed. Our work provides a complete un- 
derstanding of the structure evolution of Si on Ag(lll), 
which is desirable for fabrication of high quality silicene 
and exploring its novel physics and applications. 

Experiments and Methods Experiments were car- 
ried out in a home-built low temperature STM with a 
base pressure of 5xl0~^^ Torr. Single crystal Ag(lll) 
sample was cleaned by cycles of argon ion sputtering and 
annealing. Silicon was evaporated from a heated wafer 
(«1200 K) onto the pre-heated substrate. The deposi- 
tion rate of silicon was kept at 0.08-0.1 ML/min (here 1 
monolayer refers to the atomic density of a ideal silicene 
sheet). The STS data were acquired using a lock- in am- 
plifier by applying a small sinusoidal modulation to the 
tip bias voltage (typically 10 mV at 676 Hz). All our 
STM experiments were carried out at 77 K. 



First-principles calculations were performed within the 
framework of density functional theory (DFT) using Pro- 
jected Augmented Wave (PAW) [1^, [1^1 pseudopoten- 
tials and the Perdew-Burkc-Ernzerholf (PBE) [17 1 form 
for exchange-correlation functional, as implemented in 
Vienna ab initio simulation package (VASP) [18]. Dur- 
ing calculations, the structures were relaxed without any 
symmetry constraints using a plane-wave energy cutoff 
of 250 eV. The convergence of energy is set to 1.0x10^"* 
eV. The relaxation process continues until forces are be- 
low 0.01 eV/A. 

Results and Discussions Silicon atoms deposited on 
Ag(lll) tend to form clusters or other disordered struc- 
tures when the substrate temperature is below 400 K 
during growth (data not shown here). As substrate tem- 
perature increases to 420 K, two ordered phases form, as 
shown in Fig. 1. The less ordered phase consists of close 
packed protrusions (labeled T), and the highly ordered 
phase exhibits honeycomb structure (labeled H). The two 
phases can coexist on the surface within a large coverage 
range, from 0.5 ML [Fig. 1(a)] to 0.9 ML [Fig. 1(b)]. 

The coexistence of phase H and T indicates that these 
two phases have quite similar formation energy and sta- 
bility. One can therefore expect similarity and relations 
between the atomic structure of these two phases. In- 
deed, the high resolution STM images in Fig. 1(c) and 
(d) show that every big bright protrusion in both phase 
T and H is indeed composed of three smaller spots that 
we refer to as a "trimer" , although in phase H the trimers 
are perfectly ordered while in phase T they exhibit some 
irregularity and distortion when looked closely. The two 
phases share the same periodicity of 1.18 nm, therefore 
the density of trimers in phase H is just twice of that in 
phase T. We further performed scanning tunneling spec- 
troscopy (STS) measurements for the two phases [Fig. 
1(e)]. The dl/dV curves exhibit very similar features of 
electronic density of states (DOS), which strongly implies 
that the two phases may share some common building 
blocks in their atomic structures. In fact, in Fig. 1(c), 
there is a noticeable point showing a corner hole and six 
protrusions surrounding it - a characteristic signature of 
formation of phase H. Moreover, we notice that phase T 
prefers to form at lower Si coverage, and slightly lower 
temperature as compared with phase H. With the in- 
crease of substrate temperature and coverage, phase T 
decreases in percentage and eventually disappear com- 
pletely at 460 K. Meanwhile, phase H can spread over 
the surface if the coverage is sufficiently high. This means 
that phase H is more stable than phase T, and phase T 
can be regarded as a precursor state of phase H. 

We now face a direct question whether these two 
phases, especially the well-ordered H phase, are the theo- 
retically proposed silicene. The phase T can be excluded 
first due to the significantly lower density of Si than that 
of phase H. We notice that the periodicity of 1.18 nm, is 
almost exactly four times the lattice constant of Ag(lll)- 




FIG. 1: (color online) (a) STM image (Vup = 1.2V) of 0.5 
ML silicon atoms deposited on Ag(lll) surface at substrate 
temperature of 420 K. The areas with phase T were marked 
by "T", while the areas with phase H (with different rotation 
angles) were marked by H1-H3, respectively, (b) STM image 
(Vtip = -1.5V) of 0.9 ML silicon atoms deposited on Ag(lll) 
surface at substrate temperature of 440 K. The areas of phase 
T and H are labelled. Notably H1-H3 mark areas with phase 
H in different orientations, (c) and (d) High resolution STM 
images (8.5x8.5 nm'^, Yup = -l.OV) showing the atomic struc- 
ture of phase T and H, respectively. The blue rhombuses in 
(c) and (d) indicate the unit cells of two phases, (e) dl/dV 
spectra taken at areas of phase T (red) and H (black) respec- 
tively. The curves have been shifted vertically for clarity. 



1x1 surface, 0.29 nm, or three times the lattice constant 
of silicene, 0.38 nm. Therefore both phases H and T 
can be written as 4x4 reconstruction with respect to the 
1 X 1- Ag substrate, or 3 x 3 reconstruction with respect to 
1x1 sihcene lattice (in this paper we refer to as 3x3). 
In fact, Ag/Si system is known as a typical "magic mis- 
matched" system, such that three times the lattice con- 
stant of Si equals exactly four times the lattice constant 
of Ag. If one assume the observed H phase to be the 
theoretically proposed silicene, it is possible to obtain a 
3x3 superstructure by placing the silicene lattice in par- 



allel with the Ixl-Ag lattice. However, the crucial point 
in such moire pattern models is that: the periodicity of 
the superstructure, or essentially the moire pattern, is 
strictly linked with the relative orientations of the two 
overlapping lattices. If one obtain a 3 x 3 moire pattern in 
one orientation, it will be impossible to observe the same 
pattern in another inequivalent crystallographic orienta- 
tion. However, as we show in Fig. 1(a), we have observed 
the formation of 3x3 domains on the same Ag terrace, 
with different orientations which are obviously inequiv- 
alent. This simple experimental fact excludes the possi- 
bility that the 3x3 structure is a silicene lattice placed 
on Ixl-Ag. In another word, the 3x3 reconstruction 
should come from the structure of the overlayer itself, in- 
stead of from the commensuration between the overlayer 
and the substrate. It is, however, noticeable that the 3x3 
reconstruction is most clearly resolved in one major crys- 
tallographic orientation, while in other orientations some 
irregular distortion of the lattice is seen, which should 
come from the influence of substrate Ag lattice. 

As noted above, the periodicity of 1.18 nm is three 
times the lattice constant of Si(lll), 0.38 nm, which is 
also close to the calculated lattice constant of silicene. 
[8| Based on STM observation of the characteristic cor- 
ner hole structure, we propose a model of phase H as 
shown in Fig. 2(a). In this model, the corner holes are 
due to missing of a hexagonal silicon rings in each 3x3 
cell of a complete honeycomb silicene structure. This 
model has been confirmed by first-principles calculations. 
In the calculation the structure was modelled with low- 
buckled silicene lattice [8[ with missing silicon rings at the 
corners. The six Si atoms around the corner are hydro- 
genated. After relaxation, there are no in-plane changes 
of the position of silicon atoms, but the atoms close to 
the corner holes (red atoms in Fig. 2(a)) move upward, 
corresponding well with the trimer feature observed by 
STM. 

Based on the atomic structure of phase H, the under- 
standing of the atomic structure of phase T becomes 
straightforward. Because the STM observation shows 
that the density of Si trimers in phase T is half of that 
in phase H, we construct the model of phase T by re- 
moving half of the silicon rings in phase H, leaving only 
one hexagonal silicon ring per 3x3 unit cell, as shown 
in Fig. 2(b). This model has also been validated by 
first-principles calculations. Similar to the calculation of 
phase H, we chose low-buckled silicene rings, saturated 
by hydrogen atoms as the original structure. The cal- 
culation results show that this model is stable. Each 
trimer corresponds to a buckled silicon ring with three 
Si atoms moving upward. Such a structure can be con- 
sidered as self-assembly of hexagonal silicon rings stabi- 
lized by weak van der Waals force. Compared with the 
honeycomb arrangement of Si rings connected by cova- 
lent bonds in phase H, the weak connection of Si rings 
in phase T might explain the observed more disordered 




FIG. 2: (color online) (a) and (b) High resolution STM im- 
ages superposed with calculated model of phase H and T, 
respectively. The red and grey balls in the models represent 
buckled and unbuckled silicon atoms, respectively. The blue 
rhombus in (a) and (b)indicate unit cells as shown in Fig. 
1(c) and (d). The double arrows indicate that the rings can 
rotate randomly along their centers. 



trimer structure, as compared with the highly ordered 
trimer structure in phase H. 

When the substrate temperature during silicon growth 
reaches 480 K, the silicon structure exhibits another 
phase with obvious moire pattern, which is long range 
ordered and can spread over the whole surface as shown 
in Fig. 3(a). The orientation of moire pattern is along 
the (110) direction of Ag(lll), and the period is about 
3.8 nm. The high resolution STM image in Fig. 3 indi- 
cates that a few complete honeycomb rings with lattice 
period about 1.0 nm is observed at the bright part of 
the moire pattern, and the other parts are rather defec- 
tive and disordered. Additionally the angle between the 
direction of moire pattern and honeycomb structure is 
about 30°. The dl/dV spectra measured on this struc- 
ture shows a peak at 0.3V and a shoulder at 0.9V, which 
is distinct from that of phase T and H, and indicating an 
essentially different structure formed. 

Considering that the complete honeycomb structure 
with 1.0 nm periodicity is only observed at special posi- 
tions on surface, we proposed that this honeycomb super- 
structure consists of fragments of single layer of silicene 
with strong interaction with the Ag(lll) substrate. In 
order to clarity our supposition, first-principle calcula- 
tions have been performed. The structure model is con- 
structed as single layer, low-buckled silicene being in reg- 
istry with five Ag(lll) planes. Except for the two bottom 
Ag layers, all atoms are relaxed during the geometry op- 
timization. The energetically stable structure is shown 
in Fig. 3(d). From the calculation results we find that 
silicon atoms directly above a silver atom (red balls in 
Fig. 3(d)) are higher than other silicon atoms. As a re- 
sult these atoms should be observed as bright protrusions 
in STM image and forming a v7 x v7 superstructure 
with respect to silicene or (2-v/3x2-v/3)R30° superstruc- 
ture with respect to Ag(lll). This gives a larger honey- 



comb lattice with period 0.386 x-\/7 = 1.02 nm, in accor- 
dance well with our experimental data. The simulated 
STM image according to the calculated model is shown 
in Fig. 3(f). The similar structure features and lattice 
period as observed in experimental STM images (Fig. 
3(e)) strongly support our suggested model. Actually 
there is a slight deviation between the lattice constant of 
\/7 X \fl (1.02 nm) superstructure of silicene from that 
of 2^3x273 (1.00 nm) lattice of Ag(lll), which resuh 
in the formation of the moire pattern. The optimized 
structural model in Fig. 3(d) shows the hexagonal rings 
of silicene are twisted due to the strong interaction be- 
tween silicon atoms and silver substrate. The bright parts 
of moire pattern is where the positions of atoms in sil- 
icene are little deviated from that of Ag(lll), which make 
the honeycomb superstructure stable enough to keep the 
hexagonal rings complete. In other parts of moire pat- 
tern, the larger deviation of position between atoms in sil- 
icene and those of Ag(lll) lead to unstable honeycomb 
structure and eventually breaks the hexagonal rings of 
silicene, resulting in defective and disordered structures. 
The disordered structures were not obtained in our cal- 
culations because the unit cell we choose is much smaller 
than that of a moire pattern. The angle between the lat- 
tice direction of the \/l x \/l superstructure and (110) 
direction of Ag(lll) is 30°, so the angle between the di- 
rection of moire pattern and (110) direction of Ag(lll) 
should be zero, which has been confirmed by our experi- 
ments. 

As the substrate temperature reaches 500 K and the 
coverage is up to 0.8 ML, we observed dense honeycomb 
structure which we identify as silicene. The STM image 
in Fig. 4(a) shows a one-atom-thick silicene sheet across 
the step edges of the Ag(lll) surface without losing con- 
tinuity of the atomic lattice, which is similar to graphene 
grown on metal surfaces [19J. The high resolution STM 
image of Fig. 4(b) shows a honeycomb structure. How- 
ever, different from the reported 1x1 structure of sil- 
icene [ij] , the lattice period of the honeycomb structure 
we observed is about 0.64 nm, which is corresponding 
to a a/3 X v3 honeycomb superstructure with respect to 
the 1x1 silicene lattice. This superstructure can be ex- 
plained by a symmetric-buckled silicene model 20] which 
is shown in Fig. 4(c). The calculation of free standing 
silicene shows that the six silicon atoms in one hexagonal 
ring are not in plain: two atoms are buckled upward and 
one atom is buckled downward, forming ABA configura- 
tions. The upper buckled atoms are resolved by STM as 
the v3 X v3 honeycomb superstructure. Different from 
graphene epitaxially grown on metals 2l|, |22| , we did not 
observe moire patterns in our film, which may originate 
from the weak interaction between silicene and the metal 
substrate [23|. A typical dl/dV spectrum obtained at 
the silicene terrace (black curve in Fig. 4(e)) shows a 
shoulder at 0.3V and a peak at 0.9V, which is similar 
as the LDOS distribution measured on moire patterns 
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FIG. 3: (Color online) (a) A derivative STM image (200x200 
nm^, Vtip = 1.43V) of 0.9 ML silicon atoms deposited on 
Ag(lll) surface at substrate temperature of 480K. (b) High 
resolution STM image (15x15 nm^, Vtip = -l.OV) showing 
the atomic structure of moire patterns. The bright areas ex- 
hibit complete honeycomb rings with period of 1.0 nm while 
other areas are defective and disordered. The angle between 
the orientation of the hexagonal rings and the direction of 
moire patterns is 30°. (c) dl/dV spectra taken at the moire 
pattern phase, in which a peak at 0.3V and a shoulder at 0.9V 
are observed, (d) Calculated model of \/7x \/7 superstructure 
of silicene. The grey, yellow and red balls represent the silver 
atoms, lower silicon atoms and higher silicon atoms, respec- 
tively, (e) and (f) Experimental and simulated STM images 
(1.0 eV above Fermi energy) showing the similar structure 
features and unit cell of lattice. 



phase. Another remarkable feature is a small dip located 
at 0.5 eV which is corresponding to the Dirac point (DP) 
of silicene. The dip is not much obvious compared with 
that of graphene [2J, |25[ , which is probably due to the 
pronounced electronic DOS of the underlying Ag(lll) 
substrate superimposed on the dl/dV spectra. The devi- 
ation of the energy position of DP from the Fermi energy 
may stem from the charge transfer from the Ag(lll) sur- 
face to silicene. 

The assignment of the above phase as silicon gets a di- 
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FIG. 4: (Color online) (a) 3D STM image (30x30 nm^ Yup 
— l.OV) of a single layer of silicene island across a step edge of 
Ag(lll). (b) High resolution STM image (8x8 nm^, Yup = 
1.2V) of one monolayer silicene terrace showing the \/3 x -v/3 
honeycomb superstructure with the period of 0.64 nm. (c) 
Top and side view of schematic model of \/3 x \/3 superstruc- 
ture of silicene. The red, grey and green balls represent the 
upper buckled, in-plain, and lower buckled Si atoms, respec- 
tively. The \/3 X \/3 honeycomb superstructure is indicated 
by the black hexagon, (d) STM image (54x54 nm^, Vup = 
1.5V) of 1.2 ML silicon atoms deposited on Ag(lll) surface 
at substrate temperature of 500K showing second layer of sil- 
icene formed on the first layer of silicene. (e) dl/dV spectra 
taken on the first (black) and second (red) layer of silicene 
respectively, (f) High resolution STM image (10.5x10.5 nm'^, 
Vtip = 1.5V) of area as marked by the white rectangle in 
(d) showing atomic structure of the first and second layer of 
silicene simultaneously. 



rect proof by the observation of second layer of silicene at 
higher coverage, as shovirn in Fig. 4(d). The high resolu- 
tion STM image of Fig. 4(f) shows the atomic structure 
of the first layer and second layer silicene simultaneously. 
It is obviously that the second layer silicene also exhibit a 
V3 x vo honeycomb superstructure, which indicates that 
the v3 x v3 honeycomb superstructure should originate 
from free standing silicene and not influenced by Ag(lll) 



surface. This layer-stacked silicon structure is similar as 
graphite, and is a new structural phase of silicon, which 
may host many novel properties. The dl/dV spectrum on 
second layer of silicene (red curve in Fig. 4 (e)) resembles 
that on first layer. This striking similarity between the 
LDOS of monolayer and bilayer silicene can also confirm 
our predation that interactions between the monolayer 
silicene and Ag(lll) are as weak as that between the 
two silicene layers. 

Even at substrate temperature about 500 K during 
growth, silicon atoms tend to form the moire pattern 
phase if the coverage of silicon is considerably less than 
0.8 ML. This indicates that the atomic density of sil- 
icene is higher than that of the moire pattern phase, 
which justifies our structural model again. Increasing 
the substrate temperature to above 600 K, no structure 
of silicon can be observed anymore, leaving only a bare 
Ag(lll) surface. Furthermore, if the sample of silicene 
on Ag(lll) is annealed up to 600 K, silicene film will 
also disappear. The upper temperature limit that our 
silicene film can endure is considerably lower that that 
of graphene 2^, |26|, |27| . This is another evidence of the 
weak interaction between silicene and Ag(lll) substrate. 

Conclusion We have systematically investigated the 
structure evolution of silicene on Ag(lll). With the 
increase of the substrate temperature, silicon atoms 
on Ag(lll) overcome potential barriers and form some 
metastable structural phases such as self-assembled hon- 
eycomb building blocks (phase T and H) and incomplete 
silicene film (moire pattern structure). The most sta- 
ble phase, decoupled monolayer and bilayer silicene film, 
will appear eventually. This work provides methods to 
fabricate high quality silicene, which is essential to inves- 
tigate its novel properties, and brings it closer to the use 
in nanotechnology and other related areas. 
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